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i ABSTRACT 


This  report  describes  the  techniques  used  and  results  obtained 
in  a model  study  of  the  sensitivity  of  propagation  loss  to 
variations  in  bottom  loss  and  bottom  loss  to  variations  in 
the  geoacoustical  models  of  the  subbottom.  The  propagation 
loss  sensitivity  study  was  carried  out,  using  both  ray  and 
wave  propagation  models,  for  a variety  of  source-receiver 
depth  combinations  and  with  various  sound  speed  profile  types. 
The  bottom  loss  sensitivity  stud^-  was  carried  out  using  a 
model  developed  by  Applied  Research  Laboratories  (ARL),  Tiie 
University  of  Texas  at  Austin,  aijd  is  concentrated  on  addressing 
questions  of  hidden  depths,  substrate  ei'fects,  and  sound  speed 
and  density  profile  effects. 
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I.  INTRODUCTION 


Mathematical  models  are  useful  tools  for  predicting  underwater 
sound  propagation  and  ambient  noise  fields.  Both  the  sound  fields  and 
the  model  predictions  exhibit  sensitivity  to  environmental  variations 
and  so.irce/receiver  geometries.  To  the  extent  that  the  model  predictions 
exhibit  the  same  response  as  the  sound  fields  to  environmental  or 
geometrica].  variations,  they  are  useful  for  sensitivity  studies.  For 
such  sensitivity  studies,  one  model  input  parameter  is  varied  in  a 
systematic  manner  while  all  other  model  input  parameters  remain  fixed. 

The  resulting  variations  in  an  output  parameter  of  interest  then  describe 
the  sensitivity  of  the  model  predictions  to  the  input  parameter  variations 

Bottom  interaction  is  a complex  and  often  significant  feature  of 
underwater  sound  propagation  and  ambient  noise  fields.  While  significant 
progress  has  been  demonstrated  in  recent  years  in  predicting  ti.e 
waterborne  component  of  underwater  sound  fields,  the  bottom  interacting 
component  remains  as  a complex  challenge.  This  is  especially  true  at 
low  frequencies  when  sound  energy'  can  propagate  to  long  ranges  with 
significant  bottom  interaction.  Topographic  variations,  subbottom 
penetration  of  energ;,-  with  occasional  manifestations  of  focusing  or 
sediment  layer  interaction,  slope  enhancement,  and  seamount  blockage  or 
diffraction  all  serve  to  further  complicate  bottom  Interaction 
predictions . 

A complex  problem,  such  as  bottom  interaction  of  underwater  sound, 
is  often  best  approached  by  a systematic  examination  of  separate  features 
of  the  problem.  As  an  understanding  of  the  separate  features  is  developed 
the  complications  of  simultaneous  action  of  the  different  features  can  be 
treated.  The  present  study  has  separately  examined  the  sensitivity 


of  underwater  sound  propagation  loss  to  bottom  loss  variations,  and 
the  sensitivity  of  bottom  loss  to  sediment  geopiiysical  parameter 
variations . Results  of  these  studies  are  described  in  the  following 
two  ciiapters . 

The  results  reported  here  are  a summary  of  work  on  these  problems 
covering  the  period  November  1975  to  June  1976  under  Contract 
NCX)039-76-C-00Bl.  These  studies  are  still  in  progress  and  more  recent 
results  will  be  documented  in  other  reports  to  be  issued  under  this 
contract . 

Chapter  II  considers  the  problems  of  the  sensitivity  of  computed 
propagation  loss  to  variations  in  bottom  loss.  The  bulk  of  these 
calculations  was  carried  out  using  the  FACT  model,  althougii  the 
parabolic  equal  model  was  given  considerable  use.  In  both  cases  the 
variations  in  bottom  loss  were  restricted  to  changes  in  critical 
angle  of  a simple  phenomenological  bottom  loss  type. 

No  attempt  was  made  to  incorporate,  in  catalog  fashion,  sound 
speed  profiles  from  all  the  oceans  of  the  world.  Rather,  extensive 
use  was  made  of  profiles  typical  of  the  mid-Pacific,  North  Pacific, 
and  Indian  Ocean  regions.  These  regions  are  not  only  of  considerable 
practical  interest  but  also  their  profiles  typify  classes  Including  a 
wide  variety  of  physical  effects.  In  particular,  these  profiles  deal 
with  the  distinction  between  severe  bottom  limited  cases  and  situations 
where  propagation  with  little  bottom  interaction  is  possible. 

In  addition  to  studies  Involving  sensitivity  to  bottom  loss 
(critical  angle)  changes,  consideration  is  also  given  to  the  effects  of 
differing  water  depths  and  various  source  depths. 

Although  it  would  be  possible  to  examine  changes  in  propagation 
loss  as  a function  of  depth  at  fixed  discrete  ranges,  the  presence  of 


converg>;nce  and  shadow  zones  makes  such  an  approach  complex  and  offers 
little  practical  insight  into  the  systematic  effects  of  real  concern. 

To  eliminate  the  problem  of  dealing  with  the  full  two-dimensional 
field  (cylindrical  symmetry  is  assumed),  a range  averaged  propagation 
loss  is  used.  This  average,  described  in  detail  in  the  text,  tends  to 
smooth  the  lateral  variations  and  weight  the  influence  of  shadow  zone 
and  convergence  zone  regions  in  a useful  wai' . 

The  most  important  results  obtained  in  this  propagation  loss 
sensitivit>  study  are  itemized  below. 

1.  In  the  domain  between  the  source  depth  and  source  conjugate 
depth,  the  range  averaged  propagation  loss  displays  a regular  and 
uniform  dependence  on  critical  angle  changes. 

2.  Within  this  regular  domain  the  amount  of  bottom  influence 
on  average  propagation  loss  is  a function  of  source  and  bottom  depths. 

:5 . Outside  this  domain  of  regularity,  as  the  receiver  approaches 
either  the  surface  or  the  bottom,  the  influence  of  the  bottom  increases. 

4.  In  the  near  bottom  region,  extending  to  400  to  700  m above 
the  bottom,  the  average  propagation  loss  displays  significant  fine 
structure . 

The  second  phase  of  this  study  concerns  the  problem  of  determining 
the  sensitivity  of  bottom  interaction  to  variations  in  the  maxeup  of 
the  subbottom.  There  are  several  possible  approaches  to  this  problem, 
that  is,  it  is  possible  to  construct  various  measures  of  the  sensitivity. 
The  most  obvious  of  these  measures  includes  using  the  computed  bottom 
loss  itself,  examining  changes  in  normal  modes  and  eigenvalues,  and  using 
the  propagation  loss  in  a way  that  would  combine  the  two  aspects  of  this 
sensitivity  study. 

The  results  reported  here  are  based  on  the  first  of  these 
possibilities,  use  of  computed  bottom  loss  as  a measure  of  the  effects 
of  variations  in  subbottom  structure  and  composition.  This  method  was 
chosen  because  of  its  intuitive  simplicity,  the  direct  application  of 
results  to  a further  propagation  Loss  sensitivity  study,  and  the 
availability  of  a sophisticated  general  purpose  bottom  loss  model  at  ARL. 
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The  problem  of  dealing  with  a multitude  of  paz-ametere  that  arose 
in  the  propagation  loss  sensitivity  study  is  even  more  significant  in 
the  bottom  loss  stud>- . For  tills  reason  the  study  reported  here  focuses 
on  a small  number  of  specific  effects  and  processes.  In  this  way  an 
encyclopedic  approach  is  avoided  and  attention  can  be  directed  toward 
some  phj'sical  processes  of  importance. 

The  bottom  loss  sensitivity  study  is  centered  on  two  central 
issues:  (l)  determination  of  the  depth  in  the  sediment  to  which  any 

knowledge  of  subbottora  structure  is  necessary  (the  hidden  depth 
problem)  and  (2)  the  effects  of  variations  in  certain  acoustic 
parameters  in  the  region  above  the  hidden  depth.  The  first  of  these 
is  in  some  sense  the  most  Important  and  far  reaching  since  it  determines 
the  range  over  which  the  second  must  be  investigated. 

The  most  important  conclusions  obtained  in  this  bottom  loss 
sensitivity  study  are  the  following. 

1.  The  depth  in  the  sediment  column  to  which  detailed  knowledge 
is  necessar:/,  called  the  hidden  depth,  is  closely  related  to  the  ray 
turning  depth,  although  it  is  not  identical  with  it  except  in  the  high 
frequency  limit. 

a.  The  precise  definition  of  the  hidden  depth  (the  tolerance 
for  "hiddenness")  is  not  important;  the  hidden  depth  is  a natural  and 
pixysically  meaningful  quantity. 

b.  The  hidden  depth  is  strongly  influenced  by  the  sound 
speed  gradient  tlu-ough  its  (first  order)  effect  on  turning  depth  and 
its  (second  order)  effect  on  the  distance  below  the  turning  depth  to 
the  hidden  depth. 

c.  The  hidden  depth  lies  below  the  turning  depth  a distance 
which  is  frequency  dependent  according  to  a one- third  power  law. 

2.  In  the  low  grazing  angle  regime,  the  solid  aspects  of  a 
rock  substrate  are  generally  very  small  except  for  a narrow  angular 


regime  where  Stoneley  waves  can  caitse  a large  increase  in  loss  . 


3.  In  the  high  grazing  angle  regime,  rock  substrate  effects 
are  domitiated  by  shear  and  compressional  critical  angles. 

4.  For  sediments  of  the  silt-clay  types,  the  sound  speed 
gradient  above  the  hidden  depth  exerts  a strong  influence  on  bottom 
loss  over  the  low  grazing  angle  regime. 

5.  The  effects  of  sediment  type  variation  are  most  usefully 
studied  using  the  empirical  porosity  parajneterization  techniques 
developed  principally  by  Hamilton  at  Naval  Undersea  Center  (NUC). 
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II.  set;sitivity  of  propagation  loss  to  variations 

IN  BOTTOM  LOSS 


A.  IntrcductLon 


Fo.’  present  purposes  the  bottom  interne  i:ion  pi’oLIem  can  be 
defined  by  the  questions  of  what  influence  the  bottom  has  on  long 
range  acoustic  propagation  and  how  much  information  about  the 
(sub)bottom  is  required  in  order  to  adequately  predict  these  effects. 
Even  in  this  restricted  form  the  problem  thus  posed  is  a complex  and 
detailed  one  involving  the  iniluence  of  many  competing  processes  and 
a multitude  of  parameters.  This  complexity  has  led  to  the  necessity 
for  separating  the  problem  into  simpler  components,  a procedure  dictated 
more  by  pragmatic  computational  considerations  than  by  the  physics  of 
the  problem. 

Although  it  has  cert-ainly  been  appreciated  for  some  time  that  the 

bottom  can  exert  considerable  influence  on  propagation  loss  for  some 

configurations  and  little  inl'i.uence  in  others,  no  systematic  attempt 

has  been  made  previously  to  delineate  tliese  regimes.  The  beginnings 

1 

of  this  present  study  are  discussed  in  an  earlier  report, but  all 
the  results  given  here  are  new . 

The  propagation  loss  sensitivity  study  described  in  this  chapter 
was  designed  to  isolate  regions  in  the  water  column  of  extreme  bottom 
sensitivity  from  those  of  minimal  sensitivity.  This  isolation  can  be 
thought  of  as  determining  what  acoustic  energ>  lias  significantly- 
interacted  with  the  bottom  and  what  the  resultant  effect  on  propagation 
loss  was.  This  separation  could  therefore  be  carried  out  using  very- 
simple  phenomenological  bottom  loss  types  and  the  deviation  from  them. 
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Throughout  the  study  reported  in  this  chapter,  variations  in  bottom 
loss  were  restricted  to  changes  in  critical  angle  oy  using  a bottom 
loss  displa^ying  zero  loss  below  the  critical  angle  and  an  essentially 
infinite  loss  above.  The  use  of  such  a bottom  loss  description 
enables  the  propagation  loss  sensitivity  study  to  be  uncoupled  from 
the  bottom  loss  sensitivity  study. 

Two  propagation  models  were  used  in  this  investigation:  FACT, 

2 ’i 

a ray  trace  model  reported  by  Spofford  and  by  Baker  and  Spofford" 

4 

and  a parabolic  equation  model  reported  by  Brock.  The  nature  of  the 

ARL  implementation  of  these  models  is  described  by  Hawker,  Foreman,  and 
b 

Focke.  FACT  defines  the  bottom  in  terns  of  a plane  wave  reflection 
coefficient  applied  at  the  water- sediment  ititerfaoe.  The  parabolic 
equation  model  simulates  a critical  angle  effect  by  introducing  a 
thin  layer  below  the  water  column  having  a sound  speed  gradient  and 
zero  attenuation.  This  lay'er,  in  turn,  overlies  a homogeneous 
attenuating  layer. 

Various  specialized  terns  used  throughout  the  chapter  are 
defined  here.  Tiie  depth  at  which  the  sound  speed  is  the  same  as  that 
at  the  source  depth  Is  referred  to  as  the  source  conjugate  depth.  The 
depth  at  which  the  sound  speed  is  the  same  as  that  at  the  surface  is 
called  the  critical  depth.  The  distance  between  the  critical  depth 
and  the  bottom  (if  the  bottom  lies  below  critical  depth)  is  called  the 
depth  excess . 

J 

The  work  re.norted  here  was  carried  out  using  three  profile  types: 
Mid-Pacific,  North  Pacific,  and  Indian  Ocean.  Since  the  intent  of  this 
study  was  to  carry  out  a Uiorougb  and  systematic  sensitivity  study  for 
several  important  profile  types  and  not  to  survey  ttie  effects  of 
prt'files  fr.  m ttie  world's  t-ceans,  tt;is  selection  of  profiles  was  deemed 
adequate . 
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B. 


The  Use  of  Range  Averaged  Intensity  in  t.  Sensitivity  Study 


A sensitivity  study  is  designed  to  measure  the  Influence  of  a 
given  input  parameter  on  the  total  output . This  study  is  directed 
toward  determining  the  effects  of  bottom  loss  variations  on  the 
total  propagation  loss.  The  energy-  contributions  from  the  waterborne 
paths  dominate  the  intensity  calculations  within  the  convergence  zones, 
limiting  the  major  bottom  influence  to  the  shadow  zones. 

The  propagation  loss  sensitivity  to  bottom  loss  will  be  range 
dependent.  This  range  dependence  is  a function  of  the  remaining  input 
parameters.  Sensitivity  measured  at  designated  ranges  would,  therefore, 
be  of  very  limited  value.  Measurements  of  the  sensitivity  within  both 
the  convergence  zones  and  the  shadow  zones  would  be  very  convenient . 
Acquiring  these  results  for  a large  number  of  model  runs  would  not  be 
practica.'  because  the  location  and  widths  of  these  zones  will  vary 
from  casv.'  to  case. 

An  average  propagation  loss  calculated  over  a large  range  interval 
(covering  several  convergence  zones)  will  provide  a measure  of  the  total 
propagation  loss  sensitivity  to  bottom  loss  for  that  range  interval. 

Tlie  range  interval  between  100  and  200  nm  was  used  when  computing  the 
averages.  Propagation  within  this  range  interval  will  not  involve 
the  direct  paths  of  the  rays.  There  is  only  a nominal  dB  increase 
in  loss  because  of  the  increase  in  range  :u:ross  the  interval;  the 
run  time  for  the  parabolic  equation  model  is  not  excessive  out  to  200  nm. 
Two  forms  of  averaging  iiave  been  considered.  The  first  form  involves 
averaging  the  propagation  loss  in  decibels.  Tliis  :nethod  represents  an 
average  measurement  of  the  propagation  loss  and  emphasizes  the  iiigher 
loss  values.  The  method  that  was  finally  used  -alcuiates  the  propagation 
Joss  for  the  average  intensity  over  the  range  interval.  This  method 
represents  a physically  meaningful  average  propagation  Joss. 


For  a given  source  depth,  the  average  propagation  loss  was 
computed  at  designated  receiver  depths.  The  averages  were  then 
fitted  with  a curve  and  presented  as  a function  of  receiver  depth. 

Two  of  these  curves  will  be  presented  (Figs.  II-4  and  II-5)  along 
with  their  corresponding  maximum,  miniirium,  and  median  propagation  loss, 
both  for  a highly  reflecting  bottom  and  for  a totally  absorbing  bottom. 

In  both  cases  the  average  propagation  loss  is  within  ^ dB  of  the  median 
propagation  loss  even  though  the  differences  between  the  maximum  and 
minimum  propagation  losses  vary  greatly  between  the  two  cases. 

Therefore,  these  averages  do  present  a reasonable  representation  of 
the  propagation  losses  within  the  range  interval. 

The  use  of  similar  average  propagation  loss  versus  depth  curves 
has  been  presented  in  the  literature.  In  a study  of  the  intensity  of 
field  near  the  bottom,  Gordon  uses  an  average  propagation  loss 
defined  by  the  random  phase  mode  sum  to  investigate  the  depth 
dependence  of  propagation  loss . 

Tappert,  ' in  a discussion  of  the  parabolic  etiuation  metliod  and  its 
uses,  calculates  the  average  propagation  loss  as  a function  of  deptii.  He 
presents  a series  of  curves  for  two  source  depths  and  various  frequencies. 
He  points  out  the  surface  image  interference  and  the  near  bottom  field 
anomalies . 

C . Sensitivity  in  the  Case  of  a Mid-Pacific  Profile 

1.  General  Characteristics 


Variations  in  the  bottom  loss  were  limited  to  changes  in  the 
critical  angle,  using  values  of  0°,  1°,  ‘J)°,  10°,  lp°,  and  20°.  The  loss 
was  0 dB  below  the  critical  angle  and  50  dB  above  it.  Increases  in  the 
critical  angle  will  result  in  additional  energy  and  will  indicate  the 
Importance  of  the  higher  angular  bottom  interactions. 
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Figures  Il-i,  II-2,  and  ll-'j  deinoujtratu  the  effects  of 
changes  in  the  critical  angle  on  propagation  loss  curves  as  computed 
by  FACT.  These  curves  are  for  a mid-Paciflc  profile  with  the  bottom 
set  at  critical  depth  and  a source  depth  of  l‘j2  m (500  it).  The 
receiver  depths  are  1200  m f between  soui-ce  depth  and  its  conjugate 
depth),  2p00  m (the  source  conjugate  depth),  and  5550  m (between 
source  conjugate  and  the  bottom).  The  mid- Pacific  profile  is  shown 
in  Fig.  II-A. 


Beyond  I50  rua  the  1200  in  receiver  and  the  2p00  m receiver 
lose  the  very  deep  shadow  zones.  The  bottom  x-efracting  rays  spread 
apart  with  incx’easing  range  and  beyond  IpO  nm  the..'  are  present  at  all 
ranges  for  these  two  receivers.  This  ray  spreading  also  results  in 
the  rise  of  intensity  in  the  shadow  zones  for  the  5550  m receiver. 

Caustics  formed  below  the  soiurce  conjugate  depth  result  in 
the  strong  convergence  zone  structure  of  the  25OO  m and  5550  m receivers . 
The  propagation  loss  in  these  convergence  zones  is  independent  of  the 
bottom  loss.  Regions  outside  these  convergence  zones  do,  however,  show 
influence  of  the  bottom. 

As  the  critical  angle  decreases,  the  propagation  loss  in  the 
shadow  zones  Increases . The  total  increase  in  propagation  loss  is  a 
nominal  5 <1B  for  the  1200  m receiver  and  6 to  8 dB  for  the  2500  m receiver 
and  Is  greater  than  10  dB  for  the  5550  m receiver.  This  increase  appears 
to  be  linear  with  critical  angle  beyond  I50  nm  for  the  2500  m i-eceiver. 

For  the  5350  m receiver,  this  increase  is  greatest  between  5°  arid  10°. 

This  leads  to  0°,  10°,  and  5°  grouped  together  and  10°,  15°,  and  20° 
grouped  together. 

A detailed  study  of  propagation  loss  sensitivity  using  those 
propagation  loss  curves  would  require  a step  by  step  description  of  the 
sensitivity.  The  study  could  also  require  a separation  of  the  effects 
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within  shadow  zones  and  convergetice  zones  as  weli  as  a dex’inition  of 
the  destinction  between  these  two  zones.  Instead,  the  propagation 
losses  of  the  average  intensities  within  the  100  to  200  nm  range 
intexaral  were  calculated  for  the  sensitivity  study. 

A better  understanding  of  one  curves  resulting  from  the  use 
of  the  average  values  can  be  obtained  from  Figs.  II- A and  II-I?.  These 
figures  present  characterisitics  of  the  propagation  losses  within  the 
100  to  200  nm  range  interval.  To  facilitate  correlation  of  propagation 
loss  features  with  profile  changes,  the  mid- Pacific  profile  used  for 
these  calculations  is  shown  accompany- ing  the  propagation  loss  plot. 

This  practice  will  be  followed  throughout  this  chapter.  The  peaks  in 
the  maximum  intensity  curves  occur  at  the  source  depth  and  its  conjugate 
depth.  They  result  from  very  narrow  convergence  zones  as  evidenced  by 
the  large  separation  between  the  median  and  these  intensity  maxima. 

These  narrow  convergence  zone  peaks  influence  the  mean  by  creating 
the  minima  in  the  propagation  loss  of  the  average  versus  depth  curves 
while  leaving  the  median  virtually  unchanged.  At  other  depths  in  the 
water  column,  the  propagation  losses  are  more  unifonnly  distributed 
about  the  median  value.  This  produces  a mean  value  that  is  generally 
within  1 dB  of  the  .median  value.  The  propagation  losses  for  tlie 
average  intensities  will  therefore  represent  the  genei’al  propagation 
losses  within  the  range  interval  with  an  indication  of  the  strength 
of  the  convergence  zones  at  the  source  depth  and  its  conjugate  depth. 

These  propagation  loss  versus  depth  curves  will  be  more  useful 
with  an  understanding  of  what  a change  in  the  average  value  represents 
in  terms  of  changes  in  the  propagation  Joss  versus  range  curves. 

The  average  propagation  losses  for  the  curves  in  Figs.  II-l,  II-2, 
and  1 1- 5 are  presented  in  Table  I.  The  two  curves  for  the  12CX)  m 
receiver  (Fig.  II- l)  are  represented  by  averages  whicli  differ  by  5.6  dB. 
For  the  2^00  m receiver  the  average  propagation  losses  exhibit  a linear 
relationship  with  the  critical  angle,  approximately  a 0.7  dB  increase 
for  a decrease  in  critical  angle.  For  the  5560  m receiver  the 
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TABLE  I 

AVERAGE  PROPAGATIC’J  LOSS; 

MID- PACIFIC  PROFILE 

SOURCE  DEPTH  152  m,  BOTTOtl  DEPTH  5952  m 
MODEL:  FACT 


Receiver 
Depth  (m) 

Average  propagation  loss 
for  various  critical  anrJes 

0 

(dtig) 

1 

(deg) 

5 

(deg) 

10 

(deg) 

15 

(deg) 

20 

(deg) 

1200 

93.2 

2500 

95. b 

95.4 

95-1 

92.4 

91. b 

90.9 

5550 

100.4 

99  1 

98.2 

9o.O 

94.2 

93.2 
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largest  clianges  in  the  averages  occur  at  the  lower  arigxes . The  0“  to 
1°  change  in  critical  angle  results  in  the  largest  change  in  the 
propagation  loss  per  1°. 

Figures  II-6  and  II-7  present  the  average  propagation  losses 
for  two  range  intervals,  100  to  200  nin  and  100  to  IdO  run.  The  two  sets 
of  curves  are  seen  to  have  similar  shapes  and  separation  between 
individual  curves.  The  only  major  difference  between  the  two  sets  is 
a dB  displacement,  which  results  strictly  from  the  change  in  the 
average  range.  The  averages  from  the  100  to  200  nm  range  interval, 
therefore,  are  expected  to  be  representative  for  all  ranges . 

2 . Critical  Angle  Dependence 

The  average  propagation  losses  for  a 152  m (500  ft)  source 
using  the  mid- Pacific  profile  are  shown  in  Fig.  II-D  for  various  bottom 
loss  critical  angles.  The  propagation  loss  increases  with  increasing 
critical  angle.  There  is  a local  minimum  in  the  propagation  loss  at  the 
source  depth  and  its  conjugate  depth  for  all  critical  angles.  Maximum 
propagation  losses  generally  occur  at  blie  surface  and  at  the  bottom. 

The  water  column  appears  to  be  divided  into  three  regions : 
above  the  source  depth,  between  the  source  depth  and  the  source 
conjugate  depth,  and  below  the  source  conjugate  depth.  For  all  bottom 
types  the  propagation  loss  increases  as  the  receiver  approaches  the 
surface  in  the  first  region  and  the  bottom  in  the  third  region.  In 
the  second  region  the  maxLmum  loss  occurs  at  the  axis  depth  and  decreases 
as  the  receiver  approaches  either  the  source  or  its  conjugate  depth. 

The  differences  between  propagation  losses  for  the  0°  critical 
angle  and  the  losses  for  the  other  bottom  types  are  presented  in 
Fig.  II-8.  These  curves  present  the  intensity  contributions  due  to  the 
various  bottom  interactions.  As  the  critical  angle  increases,  the 
influences  of  the  bottom  also  increases.  In  both  regions  1 and  the 
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influence  of  the  various  bottom  interactions  increases  as  tne  receiver 
approaciies  either  the  bottom  or  the  surface.  Within  region  2,  the 
influence  of  the  bottom  is  regular,  i.e.,  independent  of  receiver  depth. 

3 • Source  Depth  Dependence 

Figures  II-9  through  11-12  present  the  average  propagation 
losses  computed  by  FACT  for  the  mid-Pacific  profile  with  a bottom 
depth  of  3952  m.  These  figures  plus  Fig.  II-6  present  a series  of 
curves  for  various  source  depths . For  each  source  depth,  propagation 
loss  curves  are  presented  for  critical  angles  of  0°  and  20°. 

The  propagation  loss  curves  exhibit  minima  at  the  source 
depth  and  its  conjugate  depth.  The  maximum  propagation  loss  between  these 
two  depths  occurs  at  the  channel  axis  (bOO  m).  The  differences  between 
the  curves  for  0°  and  20°  critical  angles  have  been  computed  arid  appear 
in  Fig.  11-13 . The  minimum  difference  increases  as  the  source  depth 
decreases.  The  amount  of  energy  represented  by  the  bottom  refracting 
rays  also  decreases  as  the  source  depth  decreases  away  from  the  axis . 

The  significance  of  the  bottom  interacting  energy  will  therefore 
increase  as  the  source  depth  decreases.  Figure  II-14  presents  these 
differences  at  a 600  m receiver  depth  along  witli  the  total  amount  of 
bottom  reflecting  energy.  (This  is  represented  as  the  total  energy  in 
decibels  relative  to  the  water  refracted  energy.)  The  bottom  refracting 
energy  is  based  on  an  omnidirectional  source  and  is  the  ratio  of  the 
source  angles  for  the  first  ra^-  to  intersect  the  bottom  and  the  ray 
intersecting  the  bottom  at  20°.  The  computed  differences  of  the 
averages  are  within  2 dB  of  the  calculated  bottom  refracting  energy. 

The  differences  in  the  averages  are  calculated  for  the  intensity  and 
are  expected  to  show  only  a reasonable  agreement  with  the  bottom 
refracting  energy  calculations. 


• O PROPAGATION  LOSS  DIFFERENCES 
BETWEEN  0°  AND  20°  CRITICAL 
angle  BOnOMS  FOR  AXIS  RECEIVER 

COMPUTED  BOTTOM  REFRACTING  ENERGY 


ENERGY  - dB 


ENERGY  - dB 
10  5 


(a)  MID-PACIFIC  PROFILE 
BOnOM  DEPTH:  3952  m 
MODEL:  FACT 
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FIGURE  11-14 

COMPARISON  OF  BOTTOM  REFLECTED  ENERGY  AND  AVERAGE 
INTENSITY  OF  THE  BOTTOM  REFRACTED  PATHS 
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The  propagation  loss  differences  in  Fig.  11-15  display  the 
regular  bottom  influences  in  the  depth  domain  between  the  source  and 
its  conjugate  depth.  Ray  theory  predicts  that  all  rays  will  pass 
through  the  depths  between  the  source  depth  and  its  conjugate  depth. 

Therefore,  within  this  domain  the  energy  contributions  from  the 
bottom  interactions  will  be  constant  with  receiver  depth.  Since 
bottom  interactions  do  not  lead  to  caustics,  their  contributions 
to  the  average  intensities  will  result  in  the  regular  dependence 
seen  in  the  model. 

The  one  exception  to  this  regular  dependence  was  the  55.5  m 
source,  where  the  near  surface  receivers  below  the  source  depth  did 
not  show  this  dependence . These  model  runs  were  made  with  surface 
image  interference.  When  this  interference  was  removed  for  the  55.5  ^ 
source,  the  differences  went  from  20  dB  to  8 dB  for  the  55*5  m receiver 
and  from  8.5  dB  to  11.6  dB  for  the  50  m receiver.  Within  the  rest  of  1 

this  domain,  the  differences  remained  11  dB.  The  only  other  major  effect 
of  removing  the  interference  was  that  the  differences  for  the  l8  m 
receiver  decreased  by  10  dB  for  all  source  depths. 

4.  Bottom  Depth  Dependence 

Receivers  between  the  source  and  its  conjugate  depth  exhibit 
a regular  dependence  on  the  bottom  loss.  Using  the  parabolic  equation 
(P.E.)  model,  average  propagation  losses  were  computed  for  the 
mid- Pacific  profile  with  various  bottom  depths  (Figs.  11-15  through 

II- l8)  to  examine  the  dependence  of  this  regularity  on  the  bottom  j 

depth.  The  propagation  losses  calculated  with  a 10°  critical  angle  | 

exhibit  only  a minimum  dependence  on  the  bottom  depth  to  within  400  m 
of  the  bottom.  As  the  bottom  depth  increases  from  ibk'J  ra  to  55^5  m, 

the  maximum  decrease  in  propagation  loss  at  any  receiver  depth  is  i 

less  than  1 dB  for  receivers  more  than  40C>  ra  off  the  bottom.  Minimum 
losses  occur  at  the  source  depth  and  its  conjugate  depth.  Between 
these  two  depths,  the  maximum  occurs  at  the  axis.  Outside  this 
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domain  the  propagation  loss  has  a generally  increasing  trend  as  the 
receiver  approaches  the  bottom  or  surface. 

Figure  II- 19  presents  the  differences  between  the 
averages  computed  with  the  1°  and  10°  critical  angles.  The  depths 
between  the  source  depth  and  the  source  conjugate  depth  exhibit  the 
regular  bottom  loss  dependence  for  all  bottom  depths.  With  bottom 
depths  1000  m greater  than  the  source  conjugate  depth,  an  increase 
of  as  much  as  2000  m in  the  bottom  depth  leads  to  a change  less  than 
1.5  dJ3  in  the  differences  between  the  source  depth  and  its  conjugate 
depth. 


The  receiver  depths  greater  than  the  source  conjugate  depth 
show  a decrease  in  the  bottom  loss  dependence  as  the  bottom  depth 
increases.  The  bottom  loss  dependence  increases  as  the  receiver  depth 
increases  to  within  'JOO  m of  tlie  bottom.  Within  400  m to  JOC  m,  the 
bottom  loss  dependence  as  a function  of  depth  varies  with  the  bottom 
depth. 

D.  Considerations  of  Additional  Profiles 


A sound  speed  profile  from  the  mid- Pacific  has  been  used  in  the 
initial  investigations  of  propagation  sensitivity.  The  profile  has  a 
single  sound  channel  with  an  axis  depth  of  600  m and  represents  the 
"classical"  profile.  Profiles  from  the  North  Pacific  and  from  the 
Indian  Ocean,  both  single  channel  profiles,  are  now  used  to  extend  the 
area  of  investigation.  The  North  Pacific  profile  has  a channel  axis 
at  120  m and  presents  a "narrow"  sound  channel  (critical  depth  2260  m). 
The  profile  from  the  Indian  Ocean  has  a channel  axis  of  I676  m and 
presents  a "broad"  sound  channel  (critical  depth  4800  m).  Again,  these 
profiles  are  shown  accompani-'ing  the  respective  propagation  loss  plots. 

Average  propagation  losses  computed  by  FACT  and  P.E.  are  shown  in 
Figs.  11-20  through  11-27  lor  the  North  Pacific  Profile.  The  avei-age 
propagation  losses  for  the  Indian  Ocean  appear  in  Figs.  11-28  through 


55 


VELOCITY  - m/sec  DIFFERENCE 


r 


HhilV 


PROPAGATION  LOSS 


in  c 
in 

o o 
-J  o 
c\j 

o z I 

CM  O O 

I ^ o 


HH  *=a: 

o 

^ ? 

^ & 
O 

CD  CX 
»-H  CX. 

u. 


o 

£X 

LU 


(T3 

> 

S. 

0) 


<v 

o> 


<C  03 
(X 


ARL  • UT 
AS-76-532 
KCF  . DR 
4 . )9  . 76 


- Hid3a 


NORTH  PACIFIC  PROFILE 

SOURCE  DEPTH:  50  m,  BOTTOM  DEPTH:  2260  m.  MODEL-  FACT 


VELOCITY  - m/sec  PROPAGATION  LOSS 

1470  1510  no  100 


o 

o 

o 

o 

8 

o 

o 

o 

o 

»o 

o 

IT) 

o 

»o 

CM 

CM 

J,  t- 

" *a:  <-> 

cC  U- 

^ o o 

3 “ c 

D-  o. 


CJ3  I— 

C CC 

CC  O 


ARL  . UT 
AS -76  -603 
KCF  - DR 
5-12-76 


lU  - Hid3Q 


8 

o 

8 


8 


o 

o 

o 

o 

o 

o 

m 

o 

<o 

CM 

CM 

I 

1 

LU 

a. 


o 

o 


E 

ez 

CO 

c 

LO 

00 

UJ 

LO 

O 

o 

CT» 

_l 

o 

4-H 

r— 

CM 

u. 

Z 

1 

o 

• • 

LT)  O 

o 

ce: 

3C 

C\J  1—4 

o 

o. 

1— 

1 h- 

n— 

O. 

•—  < 

<_) 

UJ 

t— 4 

t— 4 

Q 

< 

r— 

LL. 

uj  a. 

«T3 

H-4 

21 

cr:  o 

> 

O 

O 

r>  q: 

c 

h- 

O Q- 

<U 

o. 

»— 

4— • 

-*-> 

o 

ij_  QJ 

c 

z 

ao 

C3 

►— 4 

h“ 

< 

a: 

o: 

0> 

o 

LU 

C7^ 

z 

E 

> 

c 

c 

03 

VO 

Qi 

3: 

»— 

o. 

LlJ 

o 

LiJ 

O 

Of 

=> 

o 

oo 


O 

O 

»o 

<N 


ARL  • UT 
AS- 76- 529 
KCF  • DR 
4-19-76 


VELOCITY  - m/sec  PROPAGATION  LOSS 

1470  1510  1550  no  100 


o 

o 

o 

o 

o 

o 

>o 

o 

(Nl 

04 

N 


C 

oo 

o o 
—I  o 

CM 
z I 
CM  O O 

I H-i  O 

l-H  I— 

CJ  .. 

LU  C I— 

n:  Q.  03 

=5  0 > 

O Qi  S- 

•— » G.  <D 

Li_  +-) 

UJ  C 
•— * 

OC  01 
LU  05 


LO 

UJ 

— » LT) 

*— « CM 


111  - Hid30 


ARL  - UT 
AS-76-S31 
KCP  - DR 
4-19-76 


I 


LU 

Q_ 


Q 

O 


OO 

(/) 

O 


So 

I-.  <1: 

o 

O Oi 
>— I Q- 


C3 

<c 

or 


o 

o 

CVJ 

o 

o 


> 

s. 

<u 


0> 

CD 

c 

fO 

q; 


o 

oc 

Q_ 


o 

o 


< 

Q 


O 

CO 

00 

ro 


Q. 

LkJ 

Q 

2; 

O 


o 

CQ 


C\J 

LO 


oc 

h- 

o. 

LU 

Q 


O 

cn 


o 

oo 


ARL  . UT 
AS-  76-  540 
KCF  - DR 
4-19-76 


11-3^.  These  curves,  like  those  for  the  mid- Pacific,  show  minimuia 
propagation  losses  at  the  source  depth  and  the  source  conjugate 
depth  for  sources  outside  the  surface  duct.  The  maximum  loss 
between  these  two  depths  in  general  occurs  on  axis . 

Figures  11-23  and  II-2i+  differ  only  in  the  range  interval. 

Except  for  a 5 <TB  displacement,  the  averages  for  the  400  to  480  nm 
range  interval  present  the  same  curves  as  those  for  the  100  to  200  nm 
interval.  The  100  to  200  nm  averages,  therefore,  represent  the  general 
propagation  losses. 

1 . Source  Depth  Dependence 

The  differences  between  the  0°  and  20°  critical  angle  curves 
for  the  various  source  depths  appear  in  Fig.  11-35  f’oJ'  the  North  Pacific 
profile.  The  differences  for  depths  between  the  source  depth  and  its 
conjugate  depth  exhibit  the  regular  dependence  on  the  bottom.  There  is, 
however,  a slight  general  decrease  in  the  bottom  influence  as  the 
receiver  approaches  the  source  conjugate  depth.  For  receivers  outside 
this  domain,  the  differences  increase  as  the  receiver  depths  approach 
either  the  bottom  or  the  surface.  The  minimum  differences  for  depths 
between  the  source  depth  and  the  source  conjugate  depth  decrease  as  the 
source  depth  approaches  the  axis  depth. 

Figure  II-14  presents  the  total  amount  of  energy  in  the 
bottom  reflecting  patiis  as  a function  of  source  depth.  The  differences 
in  the  propagation  losses  for  the  axis  receiver  are  generally  within 
1.5  dB  of  this  computed  energy,  exhibiting  the  same  relationship  as 
that  seen  for  the  mid- Pacific  profile.  The  maximum  error  between  these 
two  calculations  occurs  when  both  the  source  and  receiver  are  on  axis. 
Caustic  corrections  will  dominate  the  intensity  calculations  for  this 
geometry  and  the  propagation  loss  computations  are  dominated  by  the 
bottom  refracting  rays. 
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2. 


Bottom  Depth  Dependence 


Considered  as  a futiction  of  bottom  depth,  the  propagation 
losses  computed  for  the  mid-Pacific  profile  showed  a minimal  change 
(l  dB)  in  the  bottom  loss  effects  between  the  source  depth  and  its 
conjugate  depth.  The  propagation  .loss  differences  between  10“  and 
1°  bottom  critical  angles  liave  been  computed  to  determine  the  changes 
in  the  bottom  loss  influence  for  changes  in  bottom  depth  in  both  the 
North  Pacific  and  the  Indian  Ocean  profiles.  These  curves  are 
presented  in  Figs.  11-56,  II- 37,  and  11-38. 

The  propagation  loss  differences  for  receivers  between  the 
source  depth  and  source  conjugate  depth  remain  independent  of  receiver 
depth  as  the  bottom  depth  varies  for  both  profiies.  For  an  l8  m 
source,  which  is  in  the  Indian  Ocean  surface  duct,  a decrease  of  1200  m 
in  the  bottom  depth  resulted  in  a 11  dB  increase  in  the  differences. 

Even  in  this  case,  the  differences  generally  remained  independent  of 
the  receiver  depth. 

Variations  of  300  m in  the  bottom  depth  result  in  less  than 
1 dB  change  in  the  propagation  loss  differences  when  the  source 
conjugate  depth  is  greater  than  700  m off  the  bottom.  These  variations, 
including  the  mid-Pacific  pror'ile,  were  larger  (tl  dB)  for  the  North 
Pacific  profile  and  smaller  for  the  Indian  Ocean  (±0.2  dB) . 

The  152  m source  conjugate  dejjth  was  only  200  m above  the 
bottom  when  a realistic  bottom  depth  of  3880  m was  used  with  the 
Indian  Ocean  profile.  The  propagation  loss  differences  between  tlic 
1“  and  10“  bottoms  Increased  by  5 dB. 

Receivers  below  the  source  conjugate  depth  present  bottom  loss 
dependencies  which  are  functions  of  the  profile.  For  all  three  profiles, 
a small  region  close  to  the  bottom  can  be  defined  which  shows  strong 
influences  from  the  bottom.  When  the  bottom  depth  is  less  than  the 
critical  depth,  this  region  appears  to  extend  up  about  400  m.  For 


bottom  depths  greater  than  the  critical  depth,  the  regions  appear  to 
extend  to  as  much  as  700  m above  the  bottom. 

E.  Conclusions 


The  propagation  loss  sensitivity  was  investigated  using  two 
different  models,  the  FACT  and  the  P.E.  Although  they  handle  the 
bottom  interactions  differently,  both  lead  to  similar  results. 

There  exists  a region  defined  by  the  source  depth  and  its 
conjugate  depth  where  the  propagation  loss  exhibits  a regular 
dependence  on  the  bottom  loss . The  amount  of  bottom  influence  in 
this  region  is  related  to  the  amount  of  bottom  refracting  rays . 

Since  both  the  source  depth  and  the  bottom  depth  determine  the  amount 
of  bottom  refracting  rays,  the  bottom  influence  becomes  a function  of 
these  two  depths.  Once  the  bottom  depth  exceeds  the  source  conjugate 
depth  by  1000  m,  however,  changes  in  the  bottom  depth  become 
insignificant . 

Outside  this  region  of  regular  dependence,  the  influence  of  the 
bottom  interactions  increases  as  the  receiver  approaches  either  the 
bottom  or  the  surface. 
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CHAPTER  III 

BOTTCW  LOSS  SENSITIVITY 

A,  Introduction 

The  previous  chapter  has  dealt  extensively  with  the  problem  of 
the  sensitivity  of  propagation  loss  to  variations  in  bottom  loss. 

The  present  chapter  deals  with  the  other  side  of  this  problem,  namely, 
the  factors  which  might  be  the  cause  of  particular  changes  in  bottom 
loss . More  specifically,  the  general  problem  is  one  of  determining 
the  changes  in  bottom  loss  arising  from  variations  in  the  geoacoustic 
description  of  the  bottom. 

This  problem  contains  many  factors,  the  most  important  of  which 
are  the  following: 

(1)  sediment  layering  configuration, 

(2)  sediment  type  within  each  layer, 

(5)  specification  of  sediment  attenuation  and  its  frequency 
de pendence, 

(4)  sediment  sound  speed  and  density  profiles, 

(5)  nature  and  importance  of  the  rock  basement, 

(6)  shear  wave  effects  within  sediment  layers  and  at  interfaces, 
and 

(7)  bulk  and  rough  surface  scattering. 

The  bottom  loss  model,  BOTLOSS,  developed  to  investigate  the  bottom 
loss  sensitivity  problem  (see  Hawker  and  Foreman  ) is  capable  of  treating 
all  of  these  factors  except  the  last  two.  The  effects  of  bulk  scattering 
within  the  subbottom  will  not  be  considered  in  this  report.  Rough 
surface  effects  are  to  be  considered  in  a later  report  issued  under 


this  contract.  In  the  interest  of  simplicity,  the  effects  of  sediment 
shear  waves  are  ignored,  although  they  could  be  incorporated  in  the 
bottom  loss  model  if  desired. 

The  remaining  five  problem  areas  still  constitute  a large  number 
of  degrees  of  freedom.  Consequently,  the  straightforward  encyclopedia 
approach  of  direct  parameter  variation  is  far  too  cumbersome  and  the 
problem  must  be  subdivided. 

One  obvious  restriction  is  to  impose  a limitation  on  the  range 
of  grazing  angles  to  be  considered.  If  there  is  any  bottom  loss 
whatsoever,  rays  encountering  the  bottom  at  high  angles  will  soon 
suffer  prohibitively  large  losses  because  of  their  large  number  of 
bottom  bounces.  This  effect  is  illustrated  in  Fig.  III-l,  which 
shows  the  results  of  five  FACT  runs  using  the  mid- Pacific  profile  and 
geometry  of  Fig.  Figure  IIl-l  shows  propagation  loss  versus 

range  for  critical  angles  of  0°,  p°,  10°,  15°>  and  20°.  However, 
unlike  Fig.  II-3)  the  present  case  employs  a bottom  loss  which  is 
always  0 at  0°,  rises  to  0.5  dB  at  the  critical  angle,  and  then 
jumps  to  50  dB  above  the  critical  angle,  where  it  remains  constant. 

The  zero  degree  critical  angle  case  Is  therefore  identical  with  the 
corresponding  curve  of  Fig,  II. 5.  [t  can  be  observed  that  the  bottom 
interacting  energy  which  fills  in  the  siiadow  zones  not  only  diminishes 
with  range  but  does  so  more  rapidly  for  the  energy  returned  from  the 
bottom  at  high  angles . 

On  the  basis  of  such  considerations,  it  was  decided  that  the  biilk 
of  the  bottom  loss  sensitivity  study  would  be  carried  out  over  the 
grazing  angle  range  Cfe50S25°,  Although  higher  angles  can  be  Lmportant 
for  certain  geometries,  this  angular  range  is  sufficiently  broad  to 
include  all  effects  important  in  long  range  propagation. 

A second  and  equally  important  restriction  on  the  general  problem 
is  a depth  limitation.  It  is  reasonable  to  suppose  that  there  is  a 
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FIGURE  m-1 

PROPAGATION  LOSS  versus  RANGE  FOR  VARIOUS  CRITICAL  ANGLES 
WITH  A MID- PACIFIC  TYPE  SOUND  SPEED  PROFILE 


depth  below  which  no  knowledge  of  the  subbottom  is  necessary.  Such 
a limitation  could  arise  from  attenuation  effects  or  from  the  upward 

refraction  of  a positive  sound  speed  gradient  within  the  sediments. 

Q 

This  type  of  depth  restriction  was  first  discussed  by  Williams, 
who  referred  to  it  as  the  hidden  depth. 

The  hidden  depth  question  is,  in  fact,  central  to  the  entire 
bottom  loss  sensitivity  problem  and  provides  the  most  important 
restriction  on  this  broad  problem.  Once  the  hidden  depth  has  been 
determined  for  a given  set  of  sound  speed  profiles,  sediment  types, 
etc.,  it  is  then  meaningful  to  inquire  about  the  sensitivity  of  bottom 
loss  to  changes  in  these  various  parameters  occurring  above  the  hidden 
depth. 

The  next  section  of  this  chapter  deals  in  detail  with  the  hidden 
depth  problem  and  its  dependence  on  sound  speed  profile,  frequency, 
and  sediment  type.  The  following  section  then  deals  with  some  of  the 
problems  which  occur  once  the  hidden  depth  laas  been  established.  These 
include  the  importance  of  the  rock  basement  and  its  solid  aspects 
(shear  wave  supporting  maturlal)  and  the  problem  of  sediment  partmieter 
variations . 

B.  The  Hidden  Depth  Problem 

Although  the  concept  of  a hidden  depth  appears  to  be  clear 
intuitively,  a more  precise  and  quantitative  definition  will  be 
required  if  specific  results  are  to  be  ottained. 

Q 

Williams  chose  to  approach  the  problem  via  mode  theoi’y  and 
defined  the  hidden  depth  as  the  depth  below  which  the  fractional  cliange 
in  the  eigenvalues,  caused  by  the  introduction  of  a perfect  reflector, 
was  less  than  a specifiea  tolerance.  This  approach  necessitates 
specifying  which  modes  are  to  be  considered  and  therefore  is  not 


independent  of  the  propagation  geoinetr:,’ . Williams  obviated  a 
particular  choice  of  modes  by  averaging  over  a chosen  set  of  modes. 

An  alternative  approach,  and  the  one  adopted  in  this  repox't,  is 
to  define  the  hidden  depth  in  terns  of  the  fractional  change  in  the 
bottom  loss  caused  by  the  introduction  of  a perfect  reflector  at  some 
depth.  When  this  change  is  less  than  a certain  tolerance,  the  hidden 
depth  has  been  reached.  Here,  again,  some  elements  of  the  propagation 
geonetjr,'  are  hidden  since  the  grazing  angle,  or  angular  interval,  must 
be  specified. 


Also,  in  addition  to  choosing  the  tolerance,  tliere  is  the 
possibility  of  modifying  the  definition  to  examine  the  difference 
between  a noxmial  sediment  column  atid  one  with  an  absorber  underneath, 
or  perhaps  the  difference  between  the  reflector  and  absorber  cases. 

The  most  extreme  contrast  is  obtained  by  comparing  situations  in 
which  a column  of  sediment  alternate.!^'  overlies  a perfec*  reflector 
and  a perfect  absorber.  This  definition  will  be  adopted  in  this 
report.  That  is,  the  hidden  depth,  Z^j,  is  defined  as  the  depth  at 
which  a reflector  and  then  an  absorber  can  be  inserted  in  the  sediment 
such  that  Zi(Z„)=|20  log  (|R|„  +.  /iRL.  , ) l-ej  where  |r|  is 

the  modulus  of  the  plane  wave  reflection  coefficient  and  e is  the 
tolerance.  Investigation  has  shown  that  the  precise  nature  of  the 
comparison  cases  is  not  of  great  importance  and  the  results  would  be 
little  changed  if  the  definition  were  somewhat  altered. 

In  the  remainder  of  this  section  we  shall  examine  the  hidden  depth 
problem  for  two  sediment  types,  a fine  clay  and  a coarse  silt.  The 
clac>'  sediment  has  a density  ratio  (relative  to  seawater)  of  .!  .206, 
a sound  speed  ratio  of  0.991,  and  a porosity  of  The  silt  sediment 

has  a density  ratio  of  1.59Il,  sourid  speed  ratio  of  O.966,  and  a 
porosity  of  6j5X.  The  frequency  was  chosen  to  be  90  Hz  and,  at  this 
frequency,  the  attenuations  wore  O.OO289  dB/m  for  cloy  and  0.012  dS/m 


for  silt.  The  sediment  properties  and  their  relationship  to  porosity 
follovs  the  work  of  iiuailton,^^  Akal,^^  and  Hamilton. 

Figui'es  III-2  tlu’ough  show  the  bottom  loss  versus  grazing 

angle  versus  layer  thickness  I'or  a clay  layer  overlying  a perfectly 
reflecting  surface.  These  bottom  .loss  curves,  as  well  as  all  others 
appearing  in  this  report,  were  computed  using  the  BOTLOSS  model. 

The  sound  speed  profile  was  assumed  to  be  linear  C(z)=C(0)+gz,  and 
the  gradients  g are  0.5  sec  1.0  sec  and  1.5  sec  ^ in  the  three 
drawings.  This  range  of  sound  speed  gradients  includes  those  most 
commonly  found  in  deep  ocean  sediments . 

All  three  figures  display  a prominent  diagonal  ridge  where 
orientation  with  respect  to  the  grazing  angle  axis  decreases  with 
increasing  gradient.  For  an.\'  given  angle,  at  depths  beyond  this 
ridge,  the  bottom  loss  is  seen  to  be  Independent  of  an^-  further 
increase  in  layer  thickness.  This  independence  of  layer  thickness  is 
an  indication  that,  at  this  angle,  the  hidden  depth  has  been  reached. 
Moreover,  the  decreasing  inclination  of  the  ridge  shows  that,  at  a 
fixed  angle,  the  hidden  depth  is  reached  more  quickly  for  greater  sound 
speed  gradients.  This  effect  is  caused  l)y  the  increased  upward 
refraction  of  the  sound  speed  profile  I'or  larger  gradients. 

Figures  III-5  through  III-7  show  the  quantity  1(0, z)  plotted 
versus  grazing  angle  0 and  layer  thickness  z for  this  same  clay  layer, 
with  the  same  three  sound  speed  gradients  as  in  Figs.  III-2  through 
III-4.  The  quantity  l(0,z)  is  defined  as  the  absolute  value  of  the 
difference  between  the  bottom  losses  in  the  case  of  a perfectly 
reflecting  plane  and  a perfect  absorber  inserted  at  a depth  z.  This 
difference  can  be  seen  to  be 

^ 1-20  * 20  108(|Rl^„^„,„)l 

= 20iiog(  I , 
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FIGURE  m-3 

BOTTOM  LOSS  versus  GRAZING  ANGLE  versus  LAYER  THICKNESS 
WITH  A SOUND  SPEED  GRADIENT  OF  1.0  sec-’ 
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FIGURE  nr-5 

INTENSITY  RATIO  FOR  REFLECTOR- ABSORBER  SUBSTRATES 
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FIGURE  nr -7 

INTENSITY  RATIO  FOR  REFLECTOR- ABSORBER  SUBSTRATES 
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which  is  siuiplj'  the  intensity  ratio  in  the  two  cases  expressed  in 
decibels.  The  fact  that  these  curves  liave  been  clipped  at  ^,0  dB  account 
for  the  flat  portion  of  the  surface  auove  the  px’ominent  "cliff." 

There  are  severaJ  important  points  to  be  noted  concerning 
Figs.  III-5  through  IH-7.  The  most  obvious  feature  is,  of  course, 
the  vexy  steep  "cliff"  below  which  A(y,z)  rapidly  approaches  0. 

Moreover,  if  the  curve  defining  the  bottom  of  the  cliff  were  drawn  on 
the  figures,  it  would  change  its  inclination  with  respect  to  the  angle 
axis  in  a fashion  analogous  to  the  ridge  in  the  preceding  three 
figures.  That  is,  the  larger  the  gi’adlent,  the  more  rapidly  a given 
small  value  of  A(9,z)  is  achieved. 

The  hidden  depth,  defined  b^  Is  directly  obtainable 

fx’om  these  figures  once  the  tolerance  g is  chosen.  Since  there  is  no 
particular  reason  to  ciioose  any  pax'ticular  value  for  t ( in  a 
specific  problem  € might  be  fixed),  it  might  seem  at  first  that 
a considerable  element  of  arnitrariness  remains  in  the  problem.  It  is 
just  at  this  stage,  however,  that  tixe  steepness  of  the  surface  in  the 
vicinity  of  the  cliff  plays  a x'olo . inspection  it  is  clear  that  a 

change  of  e from,  say,  0.1  dE  to  0.01  dB,  will  cause  very  little  change 
in  the  derived  value  of  the  hidden  depth.  In  other  words,  the  concept 
of  a hidden  depth  will  be  a useful  one  and  the  tolerance  g need  not  be 
specified  very  px’ecLsely. 

Figure  shows  the  hidden  depth,  defined  with  e-;0.1  dE, 

plotted  versus  grazing  angle  for  the  three  sound  speed  gradienss 
shown  in  the  preceding  groups  of  figures.  For  comparison  purposes. 

Fig.  III-9  shows  the  ray  turning  depth  versus  gi’azing  angle  for  these 
same  gradients . It  will  be  observed  that  the  hidden  depth  always 
lies  sevex’al  wavelengths  below  the  turning  depth.  This  close  associa- 
tion of  hidden  depth  and  turning  depth  is  not  surprising  since  it  is 
at  the  turning  depth  that  the  pressure  field  changes  from  oscillatory 
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FIGURE  m -9 

TURNING  DEPTH  versus  GRAZING  ANGLE  FOR  A CLAY 
LAYER  HAVING  VARIOUS  SOUND  SPEED  GRADIENTS 
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to  exponentiallj''  decaying  behavior.  In  the  investigation  carried 
Q 

out  by  Williams,  the  hidden  depth  was,  in  fact,  referenced  to  the 
turning  depth  and  was  again  found  to  lie  a feu  wavelengths  below 
it. 


Having  established  the  hidden  depth  in  one  particular  case,  we 
come  to  the  general  question  of  the  variation  of  the  hidden  depth 
with  changes  in  sediment  type,  with  different  sound  speed  profiles, 
with  frequency,  etc.  Although  these  questions  have  not  yet  been 
completely  resolved  and  these  lines  of  investigation  are  still  being 
pursued,  some  results  have  been  obtained. 

Figures  III- 10  tlu’ough  III- 12  show  tlie  bottom  loss  versus  grazing 
angle  versus  layer  thickness  for  a coarse  silt  layer  overlying  a perfect 
reflector.  These  curves  together  with  the  corresponding  ones  for  the 
perfect  absorber  are  combined  to  produce  the  difference  plots  shown  in 
Figs,  III- 13  through  III-i;^.  As  in  the  clay  layer  case,  these  two 
groups  of  figures  are  based  on  a linear  sound  speed  model  with 
gradients  of  O.p  sec  1.0  sec  and  1,5  sec 

The  same  general  features  are  present  in  these  curves  as  in  the 
clay  layer  case,  althougli  the  increased  ^amount  of  structure  tends  to 
obscure  some  trends . The  most  important  point  to  note  is  that,  once 
again,  the  concept  of  a hidden  depth  is  seen  to  be  a well-defined  one. 
The  very  abrupt  flattening  of  the  A surface  below  the  "cliff"  will  lead 
to  a very  narrow  range  of  hidden  depth  values  for  a broad  range  of 
values  of  e.  Figure  III-I6  shows  the  hidden  depth,  for  this  silt  layer, 
plotted  versus  grazing  angle  for  the  tliree  sound  speed  gradients.  Once 
again  the  hidden  depth  lies  within  several  wavelengths  of  the  tuning 
depth  (which  is  computed  according  to  Z,p=g  (cw/cos  9 - c^(o)J  where 
is  the  water  sound  speed  and  c^(0)  is  the  sediment  sound  speed  at 
the  top  of  tiie  layer). 
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Tlie  close  association  between  t]ie  hidden  and  turning  depths 

leads  at  once  to  the  conclusion  that  the  apparently  vei*y  complex 

dependence  of  hidden  depth  on  all  the  various  sediment  parameters 

is  in  actuality  not  an  overriding  problem.  Upon  writing  the  hidden 

depth  as  Z„=Z„+&X,  where  \ is  some  sound  wavelength  (say,  the  water 
n i 

wavelength)  and  6 is  a parajneter,  we  see  that  all  of  the  complexity 
is  contained  in  S alone . Moreover,  we  now  have  reason  to  believe 
that  6 is  of  order  unity,  that  is,  it  is  not  much  larger  than  one. 
in  general,  6 will  depend  on  all  the  sediment  parameters  as  well  as  on 
frequency . The  dominant  sound  speed  profile  dependence  should  be 
contained  in  the  tem  Z^. 


Thus  far  we  tiave  restricted  our  investigation  to  the  case  of  a 
linear  sediment  sound  speed  profile.  Although  it  is  true  that,  in  the 
uppemost  sediment  layers,  the  sound  speed  gradient  is  generally 
constant,  seismic  profiling  measurements  clearly  show  that  the 
gradient  decreases  with  depth.  To  gain  some  insight  into  how  such 
behavior  might  affect  the  hidden  depth  picture,  we  shall  examine  the 
exponential  sound  speed  profile  proposed  by  Williafns'’  and  given  by 
C(z)=C(0)[l^A^ie'^^-l)J‘^'^^ 


with  C(«^)=C(0)[1-A^J‘-^/^,  C'(Z}~C'(0)e"’^'^UtA^(e'^^-l)j'^/^, 


Tids  profile  is  a three- parameter  one 

and 

C'(0)=l/2  bA^C(o).  Figiu-e  III-IT  shows  C'(Z)  plotted  versus  depth  for 
various  values  of  C(o)/C(»),  with  C'(0)  chosen  to  be  unity.  It  will 
be  observed  that,  for  C(0)/C('»)  < 2/;^,  C'(Z)  < 1 for  all  depths. 

Figure  III-I8  shows  the  ray  turning  depth  for  the  case  of  water 
overlying  sediment  that  contains  such  an  exponential  sound  speed 
profile.  As  would  be  expected  from  Fig.  III-I7,  for  C(0)/C(®)  > 2/i, 
the  ray  turning  depth  will  always  be  larger  than  in  the  case  of  a 
linear  profile  (dotted  line).  For  a given  value  of  C(o)/c(6d),  the 
turning  depth  approaches  infinity  at  a certain  angle.  These  critical 
angles  are  denoted  by  the  asymptotes . 


The  most  Important  result  contained  in  these  figures  concerns  the 
deviation  of  Z_  from  linear  behavior  at  low  angles.  It  will  be  noted 
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FIGURE  m -17 

SOUND  SPEED  GRADIENT  versus  DEPTH  FOR  AN 
EXPONENTIAL  PROFILE  WITH  VARIOUS  VALUES  OF  C(0)  C(®) 


TURNING  DEPTH 


that  for  8 ^ 25°  and  0.8  2 C(oj/C('")  = 2/5,  tne  ii.aximum  deviation  from 
linear  beliavior  is  4o  ra.  For  this  range  oi'  C(g)/C(uo),  the  sound  speed 
gradient  at  200  iii  can  range  I'i'oii;  O.ut  to  0.5Y5-  Thus,  aitnougn  a 
considerable  variation  of  C'(Z)  can  occur,  tne  turtiing  depth  will 
range  only  over  a coiiiparatively  narrow  interval  . Consequently,  the 
hidden  depth  values  wnich  are  previously  determined  in  several  cases 
using  the  linear  profile  will  not  be  greatly  altered  by  the  introduction 
of  sucii  an  exponential  profile. 

As  a final  problem  in  this  hidden  depth  section,  we  sliall  briefly 
examine  the  question  of  the  frequency  dependence  of  the  hidden  depth. 
There  are  two  ways  in  which  frequency  enters  tiie  problem:  directly 

and  through  the  variation  of  attei.uation  with  frequency.  The  second 
of  these  is  the  subject  of  continuing  investigation  and  will  be 
discussed  in  a later  report;  we  shall  consider  here  only  the  direct 
frequency  dependence. 


An  estimate  of  the  direct  frequency  dependence  can  be  carried  out 

in  tiie  following  way.  First,  we  expand  the  linear  sound  speed  profile 

2 2 2 2 
wave  nuii.ber  /cj^  in  a power  series  in  Z-Z,^,  (Z)-k^2 

(at(Z-Z,j)bil/2(Z-Z^)^C-t-.  ..)  witii  a-(l+fiZ^)"^  and  b---2ll(l+pZ,p)~'^  and 

where  the  linea.'  profile  is  given  by  C^(z)=C^(l<-|3Z) . Hence,  in  the 

neighborhood  of  Z^  we  have  the  dtl'ferential  equation 


b ( Z - Z,j., ) 


p = 0 


(i) 


the  solution  of  which  is  given  by 


p_j_(Z-Z^)  = aAi 


. 1/^ 

^2(3u„^  \ 


(Zi-Z„) 


i 


+ bBi 


^2;iu 


1/3 


(2) 
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where  cos  9 and  (a)  and  (b)  are  constants.  Upon  introducing 

the  phase  velocity  v=c q/cos  0 and  the  gradient  g=CJ^^,  we  have 


(3) 


Now,  below  the  turning  point  we  must  have  a decreasing  pressure  field, 
so  we  take  b=0  and  obtain 


p^(Z-Z^)  = aAi 


(M 


Now  consider  two  frequencies  and  f^  and  suppose  that  p{Z-Z^), 
for  Z > Z_,  has  the  value  e at  f=f  at  Z=Z  . That  is. 


^’l^  = e = aAi 


03 


1/3 

\u^) 


(3) 


This  equation  determines  the  constant  (a).  We  now  require  that,  at  the 
second  frequency  f^,  the  pressure  p(Z-Z^;  f^)  equals  e at  some  depth 
Z'-Zg.  Furthermore,  we  suppose  that  Zj^=Z^+n\^,  and  we  obtain 


(V^T^ 


(6) 


which,  upon  the  introduction  of  Z2=Zr„+n6\2,  reduces  to 


(ri2jcc , 


) (?) 


/ V 0)  / 


which  yields  directiy  6=( f /l‘  . This  same  one-third  power  law  was 

n ^ 

obtained  by  Williams  in  his  investigations. 


To  recapitulate  these  results,  if  at  a frequency  f^,  the  pressure 
is  reduced  to  a factor  e at  a depth  (which  is  n wavelengths  below 
the  turning  depths),  then,  at  f=f2>  the  pressure  will  achieve  this  same 
value  e at  a distance  n6  wavelengths  below  the  turning  point.  The 
factor  6 is  then  given  by  6=(  f^/f^ )■*"'  \ This  rather  weak  frequency 
dependency  indicates,  for  exfimple,  that  the  previous  results  at  50  Hz 
would  be  altered  by  a factor  of  1,25  at  100  Hz. 


Having  investigated  the  question  of  the  hidden  depths  and  some  of 
its  determining  factors,  we  are  now  prepared  to  consider  other  sensitivity 
aspects  which  can  occur  above  the  hiddeti  depth. 

C . Bottom  Loss  Sensitivity  Problems  Above  the  Hidden  Depth 

This  section  deals  with  some  bottom  loss  sensitivity  problems  in 
the  region  of  sediment  above  the  hidden  depth.  Unfortunately,  although 
it  is  possible  to  ask,  and  to  answer,  quite  specific  questions  in  this 
domain,  tlie  ver..-  lai’ge  number  oi  reie.ani  variables  manes  systematic 
study  difficult.  Moreover,  it  is  to  some  extent  pliysically  unrealistic 
to  consider  independetit  variations  in  many  of  those  parameters  since 
the  pli;ysical  properties  of  real  sediments  couple  them  together. 


It  is  useful  to  consider  a subdivision  into  categories  which, 
while  not  completely  independent,  reflect  different  physical  processes 
and  acoustical  effects . The  most  useful  such  breakdown  is  in  terms  of 
the  following  problem  areas: 

(l)  variation  in  sediment  type  (density,  sound  speed,  and 
attenuation)  with  fixed  layering,  sound  speed  profile  function,  and 
substrate. 


(2)  substrate  effects,  such  as  the  Importance  of  sheaz’  waves 
in  the  substrate, 

(3)  sound  speed  profile  and  f^radient  effects,  with  other 
sediment  parameters  fixed, 

(4-)  effects  of  a continuously  variable  density  gradient, 

(5)  attenuation  effects  (singling  out  this  important  parameter 
that  is  difficult  to  measure  to  determine  the  required  precision  of 
for  which  detailed  knowledge 

(6)  frequency  dependence  of  all  the  above,  and 

(7)  sediment  shear  wave  effects,  treating  the  sediments  as  solids 
rather  than  fluids . 

Although  it  has  been  widely  appreciated  for  some  time  that  all 
these  various  effects  are  component  parts  of  bottom  reflection  loss, 
very  little  has  been  done  in  the  way  of  detailed  quantitative  and 
systematic  studies  of  them.  [n  large  measure,  this  is  due  to  the 
lack  of  adequate  models  incorporating  the  various  processes,  a 
deficiency  which  has  only  recently  been  remedied  by  the  BOTLOSS  model 
developed  by  ARL/UT  and  by  a model  developed  at  NUC  by  H.  Bucker. 

Some  work  in  the  direction  of  a sensitivity  study  was  carried  out  at 
SACLANT  by  St  range  rup,  although  it  was  restricted  to  isovelocity, 

isodensity  sediiaents . 

In  the  remainder  of  this  section  we  shall  consider  briefly  the 
first  four  of  the  topics  listed  above.  The  remaining  three  will  be 
considered  in  a subsequent  report  under  this  contract  as  will 
additional  and  more  complete  considerations  on  the  first  four. 
Investigation  of  the  problems  associated  with  the  first  foui’  categorie 
is  not  complete,  but  the  results  presented  here  do  indicate  trends 
and  establish  some  qualitative  bounds  on  the  importance  of  these 
questions . 


i. 


Variations  in  Sedijaent  Type 


Even  ■with  all  other  probleia  parameters  fixed  (layering, 

somd  speed  profile  function,  gradient,  and  substrate,  etc.)  the 

problem  of  the  influence  of  sediment  type  variations  on  bottom  loss 

in-volves,  at  a minimum,  variations  in  sound  speed,  density,  and 

attenuation.  One  possibility  for  circumventing  this  difficulty  is  to 

relate  these  three  parameters  to  a single  parameter  such  as  porosity 

11  12 

or  mean  grain  size.  The  work  of  Akal  and  Hamilton  has  shown 
that,  for  a considerable  variety  of  sediments,  these  three  parameters 
can  be  usefully  related  to  porosity.  Furthermore,  appropriate 
equations  have  been  developed  by  these  authors . 

Upon  fixing  all  other  paraineters  in  the  problem,  then  it  is 

possible  to  study  sediment  type  variations  by  examining  the  behavior 

of  bottom  loss  versus  angle  versus  porosity.  An  empirical  approach 

lb 

to  this  type  of  study  was  taken  by  Hall  and  Watson  who  used 
empirical  fits  to  AMOS  data  and  developed  curves  giving  bottom  loss 
versus  grazing  angle  parameterized  by  porosity. 

It  is  quite  straightforward  to  employ  the  aforementioned 
results  of  Hamilton  and  Akal  to  generate  appropriate  parameter  sets 
for  use  in  our  bottom  loss  model.  This  has  been  done  and  Figs.  III-I9 
find  III-20  show  two  resulting  three-dimensional  surfaces  in  the 
case  of  a linear  sound  speed  profile  with  a gradient  of  1.2  sec  ^ in 
a 100  m layer  and  a 200  m layer,  respectively'-.  The  substrate  is  a 
homogeneous  solid. 

There  are  two  striking  features  of  these  surfaces:  the 

bladed  structures  at  low  angles  and  the  curving  ridge  at  higher  angles 
Interpretation  of  these  curves  is  not  yet  complete,  but  it  is  certain 
that  the  prominent  curving  ridge  is  associated  with  Stoneiey  waves  at 
the  sediment- substrate  interface,  whereas  the  bladed  structures  seem 
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FREQUENCY:  50  Hi 

LAYER  THICKNESS:  200  m 

SOUND  SPEED  GRADIENT:  1.2  sec"' 


FIGURE  m-20 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
versus  POROSITY  FOR  A 200  m LAYER 
OVERLYING  A SOLID  SUBSTRATE 
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to  be  a residue  of  an  intromission  angle  efi'ect.  It  is  uotewortiiy 
that,  since  the  bladed  structures  lie  at  angles  such  tliat  the  substrate 
is  below  the  hidden  depth,  no  change  is  seen  in  these  structures  between 
the  100  m and  200  m cases . 

These  curves  show  directly  how  the  bottom  loss  can  be 
expected  to  change  with  sediment  type  variations.  Future  work  will 
include  producing  additional  such  curves  with  different  sound  speed 
gradients  and  at  different  frequencies. 

2.  Substrate  Effects 


The  most  obvious  effects  to  be  expected  from  the  presence  of 
a solid  substrate  located  above  the  hidden  depth  would  be  associated 
with  the  critical  angles  of  compressional  and  sliear  waves  in  the 
substrate . Such  critical  angle  effects  are  indeed  evident  in 
Figs.  III-21  through  III-23,  which  show  the  bottom  loss  for  a clay 
sediment  having  a constant  gradient  of  1.2  sec  and  for  the  three 
frequencies  25  Hz,  50  Hz,  and  100  Hz.  The  shear  wave  critical  angle 
is  at  55-2°  and  the  compressional  wave  critical  angle  is  at  74°.  The 
dashed  curves  show  the  bottom  loss  for  the  same  set  of  parameters 
except  that  the  substrate  shear  wave  speed  has  been  set  to  0. 

Figure  III-25  shows  an  example,  at  5°,  of  the  bladed  structure  seen 
in  the  porosity  plots  discussed  in  the  previous  subsection. 

In  addition  to  the  features  just  mentioned,  these  three 
figures  display  a very  prominent  peak  in  bottom  Joss  at  low  angles, 
much  below  the  shear  wave  critical  angle.  That  they  are  associated 
not  only  with  the  substrate  but  also  with  its  solid  aspects  is 
demonstrated  by  their  disappearance  in  the  fluid  substrate  case. 

These  peaks  are  examples  of  the  peaks  comprising  the  curving  ridge 
seen  in  Fig.  111-19 . 
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A more  dramatic  presentation  of  this  effect  is  seen  in 
Figs.  III-24  and  III-25.  These  figures  shov  bottom  loss  for  a 
vai’iable  thickness  clay  layer  overlying  a perfect  reflector  and  a 
solid  substrate  (at  the  same  depth),  respectively.  The  very  letrge 
increase  in  loss  in  the  solid  substrate  case  along  the  prominent 
ridge  is  yet  another  example  of  the  effects  of  Stoneley  waves  at 
the  interface. 

A theoretical  investigation  of  this  effect  has  been  carried 

out  and  results  have  been  obtained  showing  conclusively  that  these 

large  peaks  in  bottom  loss  are  indeed  associated  with  Stoneley  waves. 

Moreover,  the  location  of  the  Stoneley  wave  peak  can  be  predicted 

using  a relatively  simple  theory  without  the  need  for  extensive 

numerical  calculation.  These  matters  as  well  as  others  were  discussed 

1*7 

in  a previous  paper  and  are  the  subject  of  a forthcoming  paper. 

Finally,  it  is  not  suggested  that  Stoneley  waves  will 
appreciably  modify  the  concept  of  a hidden  depth.  Rather,  the 
effects  of  Stoneley  waves  are  to  lead  to  a very  large  bottom  loss 
when  the  solid  substrate  is  below  the  turning  depth  and  near  the 
hidden  depth.  That  is,  even  though  the  interface  is  near  the  hidden 
depth,  the  effect  of  the  interface  may  not  be  small. 

3 . Sound  Speed  Profile  Effect 

Figures  111-26  through  III-28  show  the  bottom  loss  for  a 
150  m clay  layer  overlying  solid  (dashed  line)  and  fluid  (solid  line) 
substrates  for  three  sound  speed  gradients.  It  is  observed  that  there 
is  always  an  angle  below  which  the  two  curves  are  coincidental. 
Moreover,  this  angle  Increases  as  the  sound  speed  gradient  increases. 
This  angle  itself  is  a manifestation  of  the  hidden  depth  effect. 
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FREQUENCY:  50  Hz 
LAYER  MATERIAL:  CLAY 
SOUND  SPEED  GRADIENT:  1.0  sec"’ 
SUBSTRATE:  SOLID 


FIGURE  ni-24 

BOTTOM  LOSS  versus  GRAZING  ANGLE  arl  - ut 

versus  LAYER  THICKNESS  FOR  A CLAY  LAYER  keh*'r®fo 

OVERLYING  A SOLID  SUBSTRATE  6-28-76 
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FREQUENCY:  50  Hi 
LAYER  MATERIAL:  CLAY 
SOUND  SPEED  GRADIENT:  1.0  sec"’ 
SUBSTRATE:  PERFECT  REFLECTOR 


FIGURE  ni-25 

BOTTOM  LOSS  versus  GRAZING  ANGLE 
versus  LAYER  THICKNESS  FOR  A CLAY  LAYER 
OVERLYING  A PERFECT  REFLECTOR 


ARL  - UT 
AS- 76- 686 
KEH  - RFO 
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FIGURE  m-26 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  A 150  m CLAY 
LAYER  HAVING  A SOUND  SPEED  GRADIENT  OF  0.5  sec-' 


ID.  1538,  1542 


FIGURE  m-27 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  A 150  m CLAY 
LAYER  HAVING  A SOUND  SPEED  GRADIENT  OF  1.0  sec-l 


AJ.SO,  the  effect  of  incx’easing  gx’adient,  at  low  grazing 
angles,  is  to  reduce  the  bottom  loss  through  increased  upward 
refraction.  At  high  angles  the  gradient  has  very  little  effect. 

The  broad  peak  in  bottom  loss  at  low  angles  (=l6°,  2^°,  and  26°)  | 

i 

occurs  at  an  angle  for  which  a ray  would  just  become  tangent  to  the  i 

substrate  interface. 

The  effects  of  a variable  sound  speed  gradient  are  shown 
in  Figs.  III-29  tlirough  III- 52.  All  these  figures  are  for  the  case 
of  a thin  4.5  m silt  layer  overlying  a thick  500  m clay  layer. 

Figure  III-29  is  for  a linear  profile  (constant  gradient)  and 
Figs.  III-30  through  III-32  are  for  three  different  exponential 
profiles . At  a given  angle  the  gradient  at  a fixed  depth  increases 
with  increasing  C(  00  )/C(0). 

It  will  be  observed  that,  below  approximately  25°,  all  four 
curves  are  identical,  thus  indicating  that  changes  in  sound  speed 
gradient  have  no  effect  on  bottom  loss . In  the  raidatigle  range  ttiere 
is  a large  effect  for  C(»)/C(o)=1.26  and  a progressively  less  effect 
for  larger  values  of  txds  paratneter.  As  C(''->')/C(0)  increases,  the 
resulting  bottom  loss  curve  is  seen  to  approach  mox’e  closely  to  tne 
linear  model  case.  This  behaviox*  is  consistent  with  the  results 
given  in  the  first  section  of  this  chapter  concerning  the  exponential 
profile . 

4 . Effects  of  a Continuously  Variable  Density 

Although  the  modifications  of  the  usual  wave  equation  (which 
were  made  because  of  continuously  variable  density)  have  been  known  since 
the  work  of  Bergmann,  no  systematic  effort  has  been  made  to  determine  the 
corresponding  effect  on  bottom  loss.  Some  general  effects  due  to  a 
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continuously  variable  density  have  been  discussed  by  Tolstoy  ’ 
although  not  in  the  context  of  this  study. 
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FIGURE  m-30 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  AN  EXPONENTIAL  PROFILE 
IN  A 500  m LAYER  WITH  C(oO)  C(o)  = 1.28 


FIGURE  III  -31 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  AN  EXPONENTIAL  PROFILE 
IN  A 500  m LAYER  WITH  C(oO)/C(o)=  1.60 
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FIGURE  m-32 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  AN  EXPONENTIAL  PROFILE 


The  bottom  loss  model  (BOTLOSS)  includes  the  possibility  of 
having  a continuously  variable  density.  Therefore  the  capability  for 
investigating  this  question  in  detail  now  exists.  Some  work  has  been 
done  in  this  direction  and  an  indication  of  its  nature  will  be  given 
here . 


By  employing  an  appropriate  transformation  it  can  be  shown 
(see  Hawker  and  Foreman^)  that  the  effects  of  a continuously  variable 
density  can  be  separated  into  three  parts: 

(1)  a refractive  effect  caused  by  a modification  of  the 
sound  speed  profile, 

(2)  a change  in  pressure  amplitude  with  depth  as  p(Z),  and 

(3)  an  interface  effect. 

Tlmt  is,  the  solution  to  the  modified  wave  equation 


(c=!(z) 


2 

k 

o 


2 

cos 


p = 0 


(8) 


can  be  written  as 


p(z)  = v/p(Z)/p(Z^)  f(z)  , (c>) 

where  Z^  is  a reference  depth  and  f(z)  is  the  solution  to  the  equation 
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Furthermore,  the  density  enters  in  the  interface  continuity  conditions, 
so  an  increase  in  density  with  depth  will  result  in  a larger  value 
of  density  being  used  in  these  conditions. 

It  is  straightforward  to  show,  using  values  of  p'/p,  that 

C(z)=C(z),  Hence  the  refractive  effects  are  negligible.  Moreover, 

examination  of  the  structure  of  the  reflection  coefficient  shows  that 

generally  the  amplitude  factor  p(z)  is  also  of  little  importance. 

The  major  effect  due  to  density  changes  is  therefore  expected  to  be 

21 

the  interface  effect . The  results  of  current  studies  of  Hamilton 
will  make  available  actual  data  concerning  the  expected  values  of  density 
gradients  in  various  sediment  types  ( in  addition  to  other  important 
parameters).  It  will  thex’efore  be  appropriate  in  the  near  future  to 
carry  out  extensive  and  realistic  model  studies  on  the  effects  of 
density  gradients  and  to  finally  resolve  this  problem.  This  is 
not  to  imply  that  this  effect  is  large.  Figure  III-53  shows  the 
bottom  loss  for  a 100  m clay  layer  with  and  without  a density 
gradient.  Although  a large  change  can  occur  in  the  vicinity  of  the 
Stoneley  wave  peak,  it  is  seen  that  generally,  and  especially  at 
low  angles,  the  density  gradient  is  likely  to  be  of  small  importance. 

Work  on  this  problem  is  continuing  and  more  definite  results  will  be 
given  in  a future  report. 


FIGURE  in  -33 

BOTTOM  LOSS  versus  GRAZING  ANGLE  FOR  A 100  m CLAY 
LAYER  WITH  AND  WITHOUT  A DENSITY  GRADIENT 


D. 


Conclusions 


The  most  important  conclusions  reached  in  this  chapter  are  the 
following. 


1.  The  depth  in  the  sediment  column  for  which  detailed  knowledge 
is  necessary,  called  the  hidden  depth,  is  closely  related  to  the  ray 
turning  depths,  aJ-though  they  are  not  identical  except  in  the  high 
frequency  limit. 

a.  The  location  of  the  hidden  depth  is  strongly  influenced 
by  the  sound  speed  gradient  through  its  determination  of  the  turning 
depth.  The  gradient  also  affects  the  distance  between  the  turning 
depth  and  the  hidden  depth,  although  this  is  less  important  than  the 
effect  on  the  turning  depth. 

b.  The  gradual  decrease  in  sound  speed  gradient  with  depth 
expected  to  occur  in  realistic  sediments  does  not  exert  a strong 
influence  on  the  hidden  depth. 

c . The  distance  between  the  turning  depth  and  the  hidden 
depth  is  frequency  dependent  according  to  a one-third  power  law. 

2.  The  most  useful  tool  for  studying  sensitivity  to  sediment 
type  variation  is  the  porosity  parameterization  developed  principally 
by  Hamilton. 

3.  Substrate  shear  wave  effects  are  small  at  low  grazing  angles 
except  where  Stoneley  waves  play  a significant  role. 

4.  The  behavior  of  bottom  loss  at  low  grazing  angles  can  be 
dominated  by  a peak  in  loss  due  to  Stoneley  waves . This  will  occur 
for  sediment  thickness  up  to  several  hundred  meters  and  primarily  for 
clay- silt  type  sediments. 


5.  Preliminary  indications  are  that  density  gradients  do  not 
yield  important  effects  although  a much  wider  class  of  examples 
needs  to  be  considered. 

6.  Apart  from  considerations  involving  Stoneley  waves,  the 
effects  of  a sound  speed  gradient  above  the  hidden  depth  arise  from 
an  interplay  between  the  upward  refractive  effects  of  the  gradient 
and  the  effects  of  absorption. 


Ill 


ACKNOWLEDGMENTS 


The  authors  wish  to  express  appreciation  for  the  material 
contributions  of  Terry  Foreman,  ARL/UT,  who  assisted  in  much  of 
the  computer  programming  used  in  this  study. 

Also,  appreciation  is  extended  to  Professor  A.  0.  Williams,  Jr, 
of  Brown  University  and  ARL/UT  and  to  Dr.  Loyd  Hampton,  ARL/UT,  for 
helpful  discussion  and  advice. 


11^ 


PT'EotDI  .0  P^GE  FIL-IED 


REFERENCES 


1.  K.  E.  Hawker,  A.  L.  Anderson,  K,  C.  Focke,  and  T.  L.  Foreman, 

"Initial  Phase  of  a Study  of  Bottom  Interaction  of  Low  Frequency 
Underwater  Sound, " Applied  Research  Laboratories  Technical  Report 
No.  76-14  (ARL-TR- 76-14),  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin,  April  I976. 

2.  C,  W.  Spofford,  "The  FACT  Model,"  Vol,  I,  AESD  Maury  Center 
Report  No.  109,  November  1974. 

3.  C.  L.  Baker  and  C.  W.  Spofford,  "The  FACT  Model,"  Vol.  II,  AESD, 

Maurj'  Center  Technical  Note  No.  TN-74-04,  December  1974. 

4.  H.  K.  Brock,  "The  AESD  Parabolic  Equation  Model,"  AESD,  Maury 
Center  Technical  Note  No.  TN-75-O7,  December  1975. 

5.  K.  E.  Hawker,  T.  L.  Foreman,  and  K.  C.  Focke,  "A  Status  Report 
on  Propagation  and  Bottom  Loss  Models  in  Use  at  ARL/UT, " Applied 
Research  Laboratories  Technical  Memorandum  No.  77-1  (ARL-TM-77-l)> 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
January  1977 • 

6.  D.  Gordon,  "Theoretical  Propagation  of  Low-Frequency  Sound  in  the 
Deep  Ocean  and  its  Interaction  with  the  Bottom"  (U),  Naval 
Undersea  Center  Report  NUC  TP  536,  January  1977-  CONFIDENTIAL 

7.  F.  Tappert,  Presentation  at  The  International  Workshop  on  Low- 
Frequency  Propagation  and  Noise,  unedited  proceedings.  Woods  Hole, 
Massachusetts,  October  1974 ■ SECRET 

8.  K.  E.  Hawker  and  T.  L.  Foreman,  "A  Plane  Wave  Reflection  Coefficient 
Model  Based  on  Numerical  Integration:  Formulation,  Implementation 
and  Application, " Applied  Research  Laboratories  Technical  Report 

No.  76-23  (ARL-TR- 76- 23),  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin,  21  June  1976. 

9.  A.  0.  Williams,  "Hidden  depths:  acceptable  Ignorance  about  ocean 
bottoms,"  J.  Acoust.  Soc . Am.  52,  1175-1179  'I  May  1976). 

10.  E.  L.  Hamilton,  "Compressional-Wave  Attenuation  in  Marine  Sediments," 
Geophysics  (620)  1972. 

11.  T.  Akal,  "The  Relationship  Between  the  Physical  Properties  of 
Underwater  Sediments  that  Determine  Bottom  Reflection,"  Marine 
Geology  3^,  (251)  1972. 


115 


PFSCSDI  G Page  BLrNK-.^OT  FIU-IED 


r 


1 


12.  E.  L.  Hamilton,  "Prediction  of  Deep-Sea  Sediment  Properties: 

State  of  the  Art,"  in  Deep-Sea  Sediments,  A.  L.  Inderbitzen 
(ed.)  (Plenum  Press,  New  York,  1S<74). 

13.  A,  0.  Williams,  Jr.,  "Acoustic  reflection  from  a structured 

sea  bottom,"  J.  Acoust.  Soc . Am.  32>  (1976). 

14.  P.  Strangerup,  "A  Detailed  Study  of  Sound  Reflections  from  a 
Layered  Ocean  Bottom,"  SACLAI'JT  Technical  Report  No.  A2, 

SACLANT  ASW  Research  Centre,  La  Spezia,  Italy,  1965. 

15.  H.  R.  Hall  and  W.  H.  Watson,  "An  Empirical  Bottom  Reflection 
Loss  Expression  for  Use  in  Sonar  Range  Prediction,"  Naval 
Undersea  Warfare  Center  Report  NUWC  TNIO,  1967- 

16.  K.  E.  Hawker,  "The  Influence  of  Surface  Waves  on  Plane  Wave  Bottom 
Reflection  Loss  for  Realistic  Ocean  Sediments,"  paper  presented 

at  the  92nd  Meeting  of  the  Acoustical  Society  of  America, 

San  Diego,  California,  November  1976. 

I".  K.  E.  Hawker,  "The  Influence  of  Stoneley  Waves  on  Bottom 

Reflection  Loss,"  in  preparation  for  submission  to  J.  Acoust. 

Soc.  Am. 

18.  P.  G.  Bergmann,  "The  wave  equation  with  a variable  index  of 
refraction,"  J,  Acoust.  Soc.  Am.  329  (19^°)* 

19.  I-  Tolstoy,  "The  Theory  of  Waves  in  Stratified  Fluids  Including 
the  Effects  of  Gravity  and  Rotation,"  Rev.  Mod.  Phys.  207 

(1963). 

20.  I.  Tolstoy,  "Effects  of  Density  Stratification  on  Solid  Waves," 

J.  Geophys.  Res.  6009  (1965). 

21.  E.  L.  Hamilton,  "Acoustic  and  Related  Properties  of  the  Sea  Floor: 
Density  and  Porosity  Profiles  and  Gradients,"  Naval  Undersea 
Center  Report  NUC  TPA59  (1975). 


116 


k. 


28  February  1977 


DISTRIBUTION  LIST  FOR 
ARL-TR-77-17 

UNDER  CONTRACT  NOOO39-77-C-OOO3 
UNCLASSIFIED 


Copy  No. 

1-5  Commanding  Officer 

Naval  Electronic  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20j6o 
Attn:  Code  520 

6 Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  A.  P.  Franceschetti 

Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
National  Space  Technology  Laboratories 


Bay  St 

. Louis,  MS  59320 

7 

Attn: 

Samuel  Marshall  (Code  3^0 ) 

8 

Herbert  Eppert  (Code  560) 

9 

Thomas  Pyle  (Code  430 ) 

10 

Hugo  Bezdek  (Code  460j 

11 

Roy  Gaul  (Code  600) 

12 

Aubrey  L.  Anderson  (Code  320 ) 

13 

Commanding  Officer 

Office  of  Naval  Research 

Arlington,  VA  22217 

Attn:  J.  B.  Hersey  (Code  102-OS) 

Commander 

Naval  Undersea  Center 
Department  of  the  Navy 
San  Diego,  CA  92152 


14  - 16 

Attn : Code  40 

17 

M.  A.  Pedersen  (Code  307 ) 

18 

R.  R.  Gardner 

19 

Edvin  L.  Hamilton 

20 

Homer  P.  Bucker  (Code  409) 

21 

H.  Morris 

22 

D.  Gordon 

117 


L 


Distribution  List  for  ARL-TR-";’7-l / under  Contract  W00039-77-C-0003  (Cont'd) 


Copy  No. 


23  - 23 

26 

27 

28 


29 

30 


Director 

Naval  Research  Laboratory 
Department  of  the  Na\’y 
Washington,  DC  20375 
Attn:  Code  3l6o 

B.  G.  Hurdle 
R.  H.  Ferris 

Ronald  Dicus  (Code  8l20) 

Naval  Oceanographic  Office 
Department  of  the  Navy 
Washington,  DC  20575 
Attn : W . H . Geddes 

K.  V.  Mackenzie 


31  Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
Liaison  Office 
Arlington,  VA  22217 
Attn:  R.  S.  Winokur 


Commander 

Naval  Air  Development  Center 
Department  of  the  Navy 
Warminster,  PA  l897^ 

32  Attn:  C.  L.  Bartberger 

33  P . Haas 

Commander 

New  London  Laboratory 
Naval  Underwater  Systems  Center 
Department  of  the  Navy 
New  London,  CT  06320 

34  Attn:  F.  R.  DiNapoli 

35  S,  R.  Santanlello 

36  R.  L.  Deavenport 

37  P.  Herstein 

38  Commanding  Officer 

Naval  Coastal  Systems  Laboratory 
Panama  City,  FL  32^1 
Attn:  E.  G.  McLeroy,  Jr. 


118 


r 

'I 

Distribution 

List  for  AHL-TR- 77-17  under  Contract  N00039-77-C-0005  (Cont'd) 

Copy 

No. 

« 

39 

Chief  of  Naval  Material 
Department  of  the  Navy 
Washington,  DC  20360 
Attn:  G.  R.  Spalding  (Code  0345) 

1 

40  - 

41 

Superintendent 

Naval  Postgraduate  School 

Monterey,  CA  93940 

42 

43 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02545 
Attn:  John  Ewing 

Earl  E.  Hays 

1 

1 

1 

44 

Morris  Schulkin,  Consultant 
9325'  Orchard  Brook.  Drive 
Potomac , 10  20854 

! 

45 

Bolt  Beranek  and  Newman,  Inc, 
50  Moulton  Street 
Cambridge,  MA  02158 
Attn:  Preston  W.  Smith,  Jr. 

46 

Science  Applications,  Inc. 
1651  Old  Meadow  Road 
McLean,  VA  22101 
Attn:  John  Hanna 

4? 

Applied  Research  Laboratory 
Pennsylvania  State  University 
P.0,  Box  50 

State  College,  PA  I68OI 
Attn:  D,  C.  Stickler 

1 

48 

Underwater  Systems,  Inc. 
5121  Georgia  Avenue 
Silver  Spring,  MD  20910 
Attn:  Marvin  S.  Weinstein 

• 

• 

49 

Geophysics  Laboratory 
Marine  Science  Institute 
The  University  of  Texas 
700  Strand 

Galveston,  TX  775pO 

Attn:  J.  L.  Worzei,  Director 

119 

> 

• ^ 

! 

L.  - 

• , 

Distribution  List  for  ARL-TR-77-I7  under  Contract  N00059"77-C-0003  (Cont'd) 


I 

I 

I 


Cony  No. 


54 

55 

56 

57 

58 


59 


60 

61 


TRACOR,  Inc. 

1601  Research  Boulevard 
Rockville,  MD  20850 
Attn:  R.  J.  Urick 

National  Oceanic  and  Atmospheric  Administration 
Atlantic  Oceanographic  and  Meteorological  Laboratory 
19  Rickenbacker  Causeway 
Miami,  FL  551^9 
Attn:  Peter  A.  Rona 

Geopljys  ical  and  Polar  Research  Center 
Department  of  Geology  and  Geophysics 
The  University  of  Wisconsin 
Madison,  WI  55700 
Attn:  C.  S.  Clay 

The  Catholic  University  of  America 
6220  Michigan  Avenue,  NE 
Washington,  DC  20017 
Attn:  H.  M.  Uberall 

Lamont-Doherty  Geological  Observatory 
Palisades,  NY  IO964 
Attn:  Henry  R.  Kutschale 

John  E.  Nafe 
R.  E.  Hout.: 

George  M.  Bryan 

Brown  University 
Providence,  RI  02912- 
Attn:  A.  0.  Williams,  Jr. 

Yale  University 

Department  of  Engineering  and  Applied  Science 
New  Haven,  CT  06520 
Attn:  Franz  B.  Tuteur 

University  of  Auckland 
Physics  Department 
Auckland,  New  Zealand 
Attn:  Alick  Kibblewhite 

Kris  Tindle 


120 


Distribution  List  for  ARL-TR-77-I7  under  Contract  N00039~T7"C-0003  (Cont'd) 


Page  No. 


62 

63 

64  - 75 


76 


77  - 78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 


Defence  Scientific  Establishment 
HMNZ  Dockyear 
Devonport,  Auckland 
New  Zealand 
Attn:  Michael  Guthrie 

R.  N.  Denham 

Commanding  Officer  and  Director 
Defense  Documentation  Center 
Defense  Services  Administration 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria,  VA  22314 

Office  of  Naval  Research 
Resident  Representative 
Room  582,  Federal  Building 
Austin,  TX  78701 

Environmental  Sciences  Division,  ARL/UT 

Glen  E.  Ellis,  ARL/UT 

Karl  C.  Focke,  ARL/OT 

Terrj'  L.  Foreman,  ARL/UT 

Ruth  Gonzalez,  ARL/UT 

Loyd  D.  Hampton,  ARL/UT 

Kenneth  E.  Hawker,  ARL/UT 

Claude  W.  Horton,  ARL/UT 

Max  K.  Miller,  ARL/UT 

Stephen  P.  Mitchell,  ARL/UT 

Steven  R.  Rutherford,  ARL/UT 

Donald  J.  Shirley,  ARL/UT 

Jack  A.  Shooter,  ARL/UT 

Steven  L.  Watkins,  ARL/UT 


121 


Distribution  List  for  ARL-TR-77-I7  under  Contract  NOOO39-77-C-OOO3  (Cont'd) 


Copy  No. 

92 

Winifred  Williams,  ARL/UT 

93 

Reuben  H.  Wallace,  ARL/OT 

94 

Library,  ARL/UT 

95  -119 

Reserve,  ARL/UT 

